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ASSESSMENT OF CONVECTION, CONDUCTION, AND 
EVAPORATION IN NUCLEATE BOILING 
by Robert W. Graham and Robert C. Hendricks 
Lewis Research Center 

SUMMARY 

Various heat -transfer mechanisms including convection, transient conduction, and 
evaporation are discussed and evaluated as to their contribution to the overall nucleate - 
boiling heat flux. Recent boiling experiments that pertain to these mechanisms are 
cited. From the evaluation, a nucleate -boiling model is proposed that includes elements 
of each of the aforementioned mechanisms. 

From a comparison of the model with experimental information on water and 
methanol, no single mechanism dominates over the entire boiling heat -flux range. It is 
true, however, that evaporation becomes the chief contributor for heat fluxes that are 
greater than 20 percent of the critical value. 

Transient conduction is an important contributor to the total heat flux over a broad 
range of conditions. The cyclic removal of heat from the surface into the liquid sublayer 
by conduction appears to be an efficient process compared with steady-state free con- 
vection. The repeated destruction of the thermal layer that ensues the departure of a 
bubble enables this transient conduction process to be repeated at a bubble site. 

The overall model presented herein emphasizes that nucleate boiling should be 
thought of primarily as a series of transient heat -transfer processes going on over 
distinct regions of the surface. 


INTRODUCTION 

Nucleate boiling is an intriguing heat -transfer process that has been the subject of 
numerous papers in the past two decades. Diverse models of the mechanism and the 
means of correlating the heat -transfer data have been presented in the literature. In 
fact, comparisons of some models and their associated heat -transfer correlations reveal 
conflicting opinions concerning the mechanisms that explain the highly efficient heat 
transport of nucleate boiling. Fairly detailed summarizations of the various nucleate- 



boiling models and associated references are presented in references 1 to 3. However, 
a synopsis of the principal models appearing in the literature is presented herein for the 
convenience of discussion. In each of these models, one heat -transfer mechanism is 
considered to be dominant. 

One of the earliest explanations of the mechanism of nucleate boiling suggests that 
the growth and departure of bubbles from a heated surface induces turbulence that 
promotes heat transfer between the surface and the liquid. The highly turbulent liquid 
transports heat into the bulk of the liquid. The bubble is a passive agent that keeps the 
fluid agitated but is never directly involved in the heat -transfer process. This model 
was the basis for several heat -transfer correlations (see ref. 4 for discussion of several 
correlations) that treated the heat transfer as if it were a turbulent, forced-convection 
process. Bubble growth rates, bubble dimensions, and bubble rise velocities were 
utilized as the convection elements in a Reynolds number definition, which was, in turn, 
a part of a Nusselt-Reynolds number correlation. Limited experimental data could be 
grouped by these parameters, but no universal correlation would apply to all fluids. It 
is interesting to observe that such Nusselt-Reynolds number correlations can be reduced 
to a form in which the heat flux is related to a temperature difference raised to some 
power (ref. 4). 

A second model of nucleate boiling that was suggested pictures the bubble as a 
piston of a pump which displaces superheated liquid away from the surface and admits 
cooler liquid to the displaced region adjacent to the surface. By this cyclic process of 
mass transfer, large quantities of heat are transported. The originators of this model 
(ref. 2) claim that the high nucleate fluxes can be predicted by such an exchange 
mechanism. This model resembles the bubble agitation model of the previous para- 
graph. 

A third model of nucleate boiling places emphasis on the evaporative process 
associated with the growth of bubbles. A correlation was developed that indicated the 
boiling flux to be proportional to the total volume of vapor generated (ref. 5). Other 
explanations of the evaporative mechanism involve a description of the evaporative 
interface that accounts for the vapor production. A microlayer of superheated liquid 
underneath the bubble was suggested by the authors of references 6 and 7. Computations 
presented in references 8 and 9 showed that idealized evaporative fluxes from the liquid 
interfaces could account for the enhanced heat transfer of nucleate boiling. 

One of the objectives herein is to examine the singular roles of convection, vapor - 
liquid exchange, and evaporation in terms of recent experimental observations of the 
nucleate -boiling process. Included are heat-transfer data and visual studies of the 
boiling mechanism. A second objective is to postulate an overall model of nucleate 
boiling that includes the composite contributions of basic heat -transfer mechanisms. 
Postulated in several recent references are nucleate -boiling models that involve more 
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than one heat -transfer mechanism in the overall model (refs. 10 to 13). The model 
„ presented herein proposes time and surface area averages of the following basic heat- 
transfer mechanisms: 

(1) Transient thermal conduction through the liquid thermal layer in the vicinity of a 

nucleation site that is preparing to bear a bubble 

(2) Evaporation from a microlayer surface underneath a bubble that is attached to 

the heater surface 

(3) Turbulent -free convection that is taking place over the surface areas not 

covered by bubbles (A zone of enhanced convection occurs in the vicinity of a 

growing bubble. ) 

Presented in this report are estimates of the various heat -transfer contributions, 
included in the model, for a series of boiling data in which detailed statistical informa- 
tion on the bubble population over a surface is available. Although several gross 
assumptions are involved in the calculations, the results illustrate how the several 
mechanisms can contribute to an overall averaged, nucleate -boiling heat flux. The 
reader should view these calculated results as a means of assessing the relative con- 
tribution of each mechanism on the overall heat transfer and not as another nucleate - 
boiling correlation. 

SYMBOLS 


A area 

a coefficient of evaporation 

C constant or coefficient 

c specific heat 

D r dry area fraction 

G mass flux 

Gr Grashof number 

g unit conversion from force to mass 
K enhancement factor 

k thermal conductivity 

M molecular weight 

Ji dimensionless ratio 

N number of sites 
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n instantaneous number of bubbles 

Pr Prandtl number 

p pressure 

Q heat rate 

q heat flux, Q/A 

R universal gas constant 

r radius 

T temperature 

U velocity 

V volumetric expansion coefficient 

a thermal diffusivity 

/3 contact angle 

A differential 

6 thickness 

0 parameter that consolidates properties for free convection 

A. heat of vaporization 

p viscosity 

p density 

t time 

# percent of heater area covered by both single and multiple bubbles (see 

table V, p. 9) 

cp m average area under multiple bubble 

<p a a average area fraction under single bubble 

b } <iV 

^s av area ^ raction influence for single bubble 

Subscripts: 

av average 

b bubble or bulk 

calc calculated 

cond conduction, transient 

conv convection 
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cr 

* enh 
evap 
exp 
gr 
h 
l 

m 

prep 

r 

ref 

s 

sat 

v 

w 


critical 

enhanced 

evaporation 

experimental 

growth period 

heater 

liquid 

mean, average or merging (bubble) 

preparation 

receiver 

reference 

single (bubble) 

saturation 

vapor 

waiting period or wall 


DISCUSSION OF MECHANISMS 

In this section, the various boiling models that have been suggested as the reasons * 
for the large heat -transfer coefficients of nucleate boiling are discussed individually in 
the light of experimental evidence. Schematic diagrams of the models are shown in 
figure 1. The possibility of any one of these models being the dominant model in 
nucleate boiling is evaluated. 


Bubble Agitation Mechanism 

Considerable experimental evidence shows that an appreciable degree of fluid mixing 
occurs near the heater surface during boiling. As portrayed in figure 1(a), the presence 
of bubbles growing on the heater surface and departing therefrom induce this mixing or 
agitation process. In reference 14, schlieren and shadowgraph high-speed motion 
pictures were made of the thermal boundary layer in the vicinity of a bubble site. From 
this study, it was clear that the development of a bubble caused appreciable agitation in 
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the thermal layer. As the bubble grew, it appeared to displace the thermal layer, and 
agitation was visible within the adjacent thermal layer as far away as one bubble 
diameter from the center of the growing bubble. The surface area subject to this agita- 
tion is defined as the area of influence (see ref. 14). 

Figure 2 presents three frames of high-speed motion pictures showing the growing 
bubble and the displaced surface thermal layer that extends outside the bubble. A more 
vivid picture of this phenomenon can be obtained by viewing the film supplement to 
reference 14. This area of influence of a growing bubble was observed in a different 
manner by Gaertner and Westwater (ref. 15), who observed the density of bubble sites 

across a boiling surface by having each bubble ’’print” its presence on the surface by a 

* 

plating technique. They observed that the packing density of the bubbles produced an 
area of influence which was equivalent to that observed in reference 14. Apparently, the 
agitation destroys the thermal layer and limits the surface area that can be covered 
with bubbles. 

The real question concerning the agitation phenomena in nucleate boiling relates to 
their influence on the level of the heat-transfer coefficient. Certainly, fluid oscillation 
or translational motion produces enhancements but are these enhancements sufficient to 
raise the heat-transfer coefficient by 1 or 2 orders of magnitude beyond free convection? 
In reference 10, Kast states that the fluid friction present would prevent the development 
of convective velocities of sufficient magnitude to produce heat -transfer coefficients 
equivalent to nucleate -boiling values. His analytical arguments were based on a forced- 
convection, Nusselt correlation of the boiling phenomenon. Kast cited some work 
wherein he measured the convective effects of bubbles in a bubble tower and observed 
maximum coefficients of 2000 watts per square meter per °C, which are 2 orders of 
magnitude below that observed with nucleate boiling. The author claimed that the bubbles 
left the surface in the same manner as they did in pool boiling. 

In reference 16, the authors studied the pool heating of hydrogen in both its sub- 
critical and supercritical pressure states. They observed that, in the supercritical 
state, a boiling-like action took place at the heating surface. Agglomerates of warm 
fluid rose from the surface with the appearance of bubbly vapor production. The pool 
appeared to be highly agitated. Nevertheless, the heat -transfer coefficients for the 
supercritical fluid were almost 2 orders of magnitude less than those observed for 
nucleate boiling. The same kind of result was noted for water (ref. 17). 

Forster and Greif (ref. 2) pointed out that the bubble -agitation model would imply 
that the heat flux is strongly dependent on the temperature difference between the wall 
and bulk temperatures. Experimental evidence shows that the heat flux correlates 
better with T w - T gat for a range of subcooling. Employing T w - T b results in 
families of curves separated parametrically by the subcooling. 

From the evidence cited in this section, bubble agitation does not appear to be a 
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singular cause of the large heat -transfer coefficient noted in nucleate boiling. The 
bubble -agitation mechanism can be an important contributor to the effective nucleate - 
boiling coefficient. The enhanced contribution of bubble agitation to a free -convection 
process was studied in reference 18. Gas bubbles generated by electrolysis at a heat- 
transfer surface (no heat of vaporization involved) produced up to fourfold enhancement 
in the free -convection process. The degree of enhancement was primarily dependent 
on the wall- to bulk-temperature difference and secondarily dependent on the electrolytic 
rate of gas -bubble production. The enhancement in the heat-transfer coefficient was 
greatest at the lower temperature differences. To maintain only free convection, the 
liquid (sodium hydroxide) was highly subcooled. 


Vapor-Liquid Exchange Mechanism 

Forster and Greif (ref. 2) proposed a model of nucleate boiling known as the vapor- 
liquid exchange mechanism. This model, while similar in some respects to the bubble - 
agitation model, was devised in response to certain objections, as discussed in the 
previous section. 

The vapor -liquid exchange mechanism (fig. 1(b)) is a means of pumping a slug of 
hot liquid away from the wall and replacing it with a cooler slug. The cooler slug picks 
up heat from the wall before being transported out into the bulk. A growing and depart- 
ing bubble acts as the piston in this pumping process. In discussing the merits of this 
model, Forster and Greif contended that vaporization could not explain the high heat 
flux associated with nucleate boiling. They argued that the enthalpy -carrying capacity 
of the liquid slug was the answer and formulated a ratio which is the ratio of the 
heat capacity of a liquid slug to the latent heat content of a comparable volume of vapor. 
This ratio is expressed in equation form as 

cp, 

Jl = — — AT (1) 

AP V m 

where AT m is the difference between film temperature and bulk temperature. Because 
this ratio was computed to be of the order of 10^ from data obtained by Gunther and 
Kreith (ref. 19), Forster and Greif concluded that vaporization could not be important. 
However, serious objections can be raised as to how Forster and Greif estimated the 
vaporization contribution. The volume of the departing bubble was assumed to be a 
measure of the evaporative process, and a simultaneous evaporation -condensation 
process was not considered to take place within the growing bubble. They did admit 
such a possibility by referencing the early work of Snyder (ref. 6), which had just been 
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started when their paper was published. In the next section of this report, which deals 
with evaporative mechanisms, the work of Snyder and others is discussed. 

One of the basic assumptions of the liquid -vapor -exchange calculation supposes that 
a volume of liquid the size of a departing bubble is at an average temperature between 
the wall and the bulk values. This assumption was questioned by Bankoff (ref. 1) who, 
in particular, was critical of assuming that the hot-liquid slug was at film temperature, 
which would imply a thermal-layer thickness of the order of the maximum bubble radius. 
Measurements of the thermal-layer thickness that were made by probing (ref. 20) and 
by visual techniques (ref. 21) indicate that the bubble radius is generally much larger 
than the thermal layer. However, shadowgraphs and schlieren photographs of the 
thermal layer (ref. 14) did give credence to that aspect of the model which pictured a 
cooler slug of liquid being heated adjacent to the hot wall. The waiting period of the 
ebullition cycle occurs during this process. It is important to emphasize that the liquid 
being heated was within the thermal layer and that it was a thin layer of liquid that 
heated cyclicly during ebullition. In reference 22, Hsu pictured the development of the 
thermal layer in the vicinity of a bubble site as a prerequisite to bubble incipience. In 
his analysis, the thermal layer was treated as a semiinfinite slab undergoing transient 
conduction. 

Summarizing the discussion of the vapor-liquid exchange mechanism, it has been 
shown that there is experimental proof that such a mechanism exists physically. It is 
doubtful that its potential heat -flux contribution was correctly computed in reference 2 
for the reasons cited previously. The mechanism is probably significant in replenishing 
the thermal layer with cold liquid after a bubble departs, 


Evaporative Mechanism 

Any explanation of nucleate boiling must be capable of explaining the abrupt en- 
hancement in heat transfer observed when boiling begins. The two models discussed 
thus far (bubble agitation and vapor-liquid exchange) have proposed such an explanation 
in terms of the intensification of free -convection processes that are in effect at heat 
fluxes below the boiling incipience value . A third explanation is that the abrupt change 
in the boiling-heating curve is attributable to evaporation. 

Evaporation from a free liquid surface can be considered as an effusive molecular 
flow (ref. 23) and thus can be treated as a problem in the kinetic behavior of a gas or a 
vapor. In this section, the rate equation for evaporation is reviewed, and numerical 
estimates of heat -flux values possible from evaporating surfaces are presented. 

Consider a free liquid surface at saturation conditions and above the liquid a blanket 
of vapor of the same substance as the liquid. This vapor is also at saturation condi - 
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tions. A kinetic equilibrium exists whereby a certain rate of evaporating molecules 
from the liquid will be replaced by an equal rate of condensing molecules from the 
vapor. From kinetic theory, the mass rate of molecules involved in this evaporation- 
condensation exchange is given by 


r27rRT 


G = 


sat 


- 1/2 


g c M 


(2a) 


In a real boiling process in which a net production of vapor is involved, the afore- 
mentioned kinetic equilibrium is upset by a temperature or pressure difference, which 
lowers the pressure in the vapor with respect to the liquid pressure. If the vapor is at 
some receiver pressure p r , the net evaporation rate is given by the following equation, 
which is a modification of equation (2a): 


G 


evap 


, 2 ' R V '‘ 1/2 

g c M 


(P; - P r ) 


(2b) 


It has been observed experimentally that this rate equation is correct for atomic or 
molecular structures that exhibit geometric symmetry, but comparatively complex 
molecules like water or long-chain hydrocarbons exhibit evaporation rates that are a 
small fraction of those predicted by equation (2b). The ratio of the measured to the 
theoretical evaporation rate is defined as the evaporation coefficient a. Wyllie (ref . 24) 
measured the coefficients of evaporation for several fluids, some of which are given in 
table I. These coefficients were obtained by measuring the mass flow through a diffusion 
hole into a hard vacuum system and by comparing the experimental result with evapora- 
tion rates predicted from equation (2b) (with the 
receiver pressure equal to zero). Wyllie (ref. 24) 
envisioned a three -layer model of the evaporation 
mechanism. First was the liquid substrate with a 
compact, ordered layer of liquid molecules. 

Above that was a mobile absorbing layer of 
molecules. On top of the mobile layer was the 
vapor. For a high coefficient of evaporation, the 
rate of exchange between the liquid substrate and 
the mobile layer must be much greater than that 
from the mobile layer into the vapor. It is inter- 


TABLE I. - COEFFICIENTS OF 


EVAPORATION (REF. 24) 


Fluid 

Coefficient of 
evaporation, 
a 

Carbon tetrachloride (CC1 4 ) 

1 

Benzene (CgHg 

.90 

Chloroform (CHClg) 

.16 

Ethanol (C^OH) 

.02 

Methanol (CHgOH) 

.045 

Water (HgO) 

.040 
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esting to observe (as was pointed out in ref. 24) that those fluid molecules with large 
dipole moments had lower rates of evaporation. For example, carbon tetrachloride 
(CC1 4 ) with a zero dipole moment had an evaporative coefficient of 1, whereas water 
with a dipole moment of 1. 85 had an evaporative coefficient of only 0. 04. Other values 
of evaporative coefficients appear to correlate (inversely) with the magnitude of the 
dipole moment. Reference 25 presents a more recent compilation of evaporative coeffi- 
cients and an excellent review of the topic. Refined experiments for several fluids made 
by Eagleton and his colleagues (refs. 26 and 27) corroborate the general magnitude of 
the values published in reference 24. 

In several papers, the mass transfer and heat removal during surface evaporation 
was computed from the kinetic theory. Plesset (ref. 8) computed the net flux of vapor 
being transported from a hot, liquid interface to a cooler, liquid surface. The Plesset 
equation included a condensation coefficient as an unknown, which could be computed by 
the use of the equation. Usually, for equilibrium evaporation, the condensing coefficient 
is assumed to be equal to the evaporative; but for nonequilibrium situations, this need 
not be true. 

Bankoff (ref. 9) made use of reference 8 to show that the upper limit of nucleate 
boiling for several hydrocarbons was predictable as an evaporative-condensing process. 

It is interesting to examine the theroetical heat -transfer rate that is possible 
through evaporation of a fluid into a receiver. When the heat of vaporization is taken into 
account, equation (2b) can be modified to express the flux of heat leaving the surface. 

The coefficient of evaporation a is included in the equation: 


q = aX( 


f2ffRT 


sat 


- 1/2 


g c M 


<P; - P r ) 


TABLE II. - THEORETICAL EVAPORATIVE HEAT-FLUX RATES 
[Saturation pressure, 1 atm.] 


Fluid 

Saturation 

temper- 

ature, 

T sat> 

°R 

Heat of 
vaporization, 
A, 

Btu/lb mass 

Molecular 

weight, 

M 

Coeffi- 
cient of 
evapora- 
tion, 
a 

Evaporative heat flux 

Btu/ (hr)(ft^) 

W/m 2 

Mercury 

1135 

126 

200.6 

1 

44. 4X10 6 

140X10 6 

Water 

672 

970 

18 

.04 

2. 78 

8.8 

Ethanol 

632 

363 

46.1 

.02 

.85 

2.7 

Methanol 

608 

473 

37.0 

.045 

2.1 

6.6 

Nitrogen 

139 

86 

28.0 

1 

4.7 

14.8 

Hydrogen 

36 

190 

2 

1 

20.4 

64.5 


(2c) 
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Equation (2c) was employed in computing the rates of heat transfer for mercury, 
water, ethyl alcohol, liquid nitrogen, and liquid hydrogen. These rates are compared in 
table II at a saturation pressure of 1 atmosphere and with the vapor being drawn into a 
perfect vacuum (p f = 0). The results in the table are highly idealized and represent the 
case of a steady-state evaporation into a vacuum reservoir. This case is not represent- 
ative of a nucleate boiling process where (1) the process is not steady, (2) the receiver 
pressure is not a vacuum, and (3) evaporation does not occur over the entire heating 
surface simultaneously. Nevertheless, table II does show that evaporation is a highly 
efficient heat -transfer mechanism that can be 2 orders of magnitude greater than liquid- 
free convection. Consequently, the presence of evaporation over only a small portion 
of a heat-transfer surface at any one time could increase the overall heat -transfer rate, 
which could be the explanation of the almost discontinuous jump in heat transfer noticed 
when nucleate boiling begins. 


Surface Evaporation in Nucleate Boiling 

The evaporative mechanism has been demonstrated analytically to account for 
sizeable heat-flux rates. The problem of showing experimentally that a surface - 
evaporation phenomenon is present with nucleate boiling, remains. A number of inde- 
pendent investigators (refs. 3, 7, and 28 to 30) noted appreciable diminution of local 
wall temperatures underneath a bubble during the growth period. Figure 3 (from 
ref. 3) graphically illustrates this observation. The figure includes the transient 
surface -temperature history underneath a bubble and a computed curve of the transient 
heat flux during the ebullition cycle. High-speed motion pictures of the bubble history 
show the condition of the bubble at various times in the transient. The heat -flux curve 
points to the importance of the comparatively short interval of rapid bubble growth in 
contributing to the overall heat flux. As indicated by the temperature curve of figure 3, 
the severe drop in local wall temperature during bubble growth appears to contradict 
the liquid-vapor -exchange concept, which presumes that surface cooling would not start 
until after the bubble (or bubbles) departed and quenching had begun. It also contradicts 
a presumption of the bubble stirring, or agitation model, which would envision appreci- 
able cooling of the local surface underneath a bubble because of induced turbulence as 
cooler liquid replaces the bubble volume. 

Estimates of the heat -transfer rate during bubble growth are presented in refer- 
ence 3. With water as the fluid, at somewhat subcooled conditions, the heat -transfer 
rates were as high as 500 000 Btu per hour per square foot (1 580 000 W/m ). This 
maximum value approaches the magnitude of the idealized heat flux predicted from the 
kinetic theory of evaporation (see table II, p. 9). Mason and Bankoff (ref. 31) reported 
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values of heat-transfer coefficient at a surface where rapidly growing and collapsing bub- 
bles were present. These coefficients ranged in magnitude from 13 000 to 300 000 Btu per 
hour per square foot per °F (74 000 to 1 700 000 W/(m^)(°C)). Such a range of heat- 
transfer coefficient could produce heat-flux rates comparable to, or even an order of mag- 
nitude higher than, those reported in reference 3. Direct comparison between refer- 
ences 3 and 31 cannot be made because the fluid-state conditions were not comparable. 

The inclusion of the concept of an evaporative mechanism at the base of a bubble 
necessitates the presence of a thin liquid layer there. Such an assumption was made by 
Snyder and Edwards (ref. 6) and also by Moore and Mesler (ref. 7) in proposing the 
presence of a microlayer of liquid underneath the bubble. Experimental verification of 
the presence of the microlayer has been established by at least two independent investi- 
gations. Hospeti and Mesler (ref. 32) used a chemical deposition technique involving 
calcium sulfate which contained radioactive sulfur 35. The radioactive sulfur made 
possible the measurement of the amount of deposit residue left by a sequential history 
of many bubbles developing at a particular site. Autoradiographs of the deposits could 
be related to the amount of liquid boiled off in generating the bubbles; thus, estimates 
of the thickness of the liquid microlayer under the bubble could be made. For average 
heat fluxes ranging from 7000 to 20 000 Btu per hour per square foot (22 000 to 63 000 

o r* A 

W/m ), the microlayer thickness was observed to vary from 19x10 to 10 inch 
(4.84X10 -7 to 2. 54x10 ~ 6 m). 

Sharp (ref. 33) employed three kinds of experimental techniques to detect the presence 
and magnitude of the microlayer. * The simpler of the techniques was the use of an 
electrical continuity probe that was inserted into the base of a growing bubble. The indi- 
cation of electrical continuity between the probe and the heating surface was interpreted 
as the presence of a thin liquid layer. The more sophisticated method involved the 
application of the optical principles of interference fringes (Newton’s rings) to measure 
the thin microlayer. Sharp was able to make estimates of the local thickness of the 
microlayer at various radial positions beneath the bubble at several time intervals of 
bubble growth. Curves of microlayer thickness from reference 33 are reproduced in 
figure 4. It is apparent that the measured microlayer thickness varied between 1000 
and 4000 A, or approximately 4xl0“ 6 to 16xl0 -6 inch (1. 02xl0 -7 to 4. 06xl0 -7 m). These 
values are somewhat smaller than those observed by Hospeti and Mesler (ref. 32), but 
it is remarkable that there is an overlapping range of general agreement. Sharp used 
the measured values of the depletion rate of the microlayer to predict an average heat- 
flux rate underneath the bubble. The average rates varied between 30 000 and 50 000 
Btu per hour per square foot (95 000 to 158 000 W /re?) for the area under the bubble. 

*A NASA-Lewis motion picture (C-252) describing the experimental methods and in- 
cluding film clips from the research sequences is available on loan. A request card and 
description of the film are included at the back of the report. 
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PROPOSED MODEL OF NUCLEATE- BOILING MECHANISM 


As a result of the evaluation of a number of mechanisms of nucleate boiling pre- 
sented in the previous sections, the following overall model of the nucleate boiling cycle 
was constructed. This proposal incorporates specific features from the several models 
that appear to agree with recent experimental evidence. 


Boiling Surface 

Bubble area . - A boiling surface such as that depicted in figure 5(a) contains N 
number of boiling sites that are active at a given heat flux. Instantaneously, there are 
bubbles, which can be single, or multiple, over a certain fraction of these sites. 

For an intermediate range of nucleate-boiling heat flux, both single bubbles and multiple 
bubbles that consist of several single bubbles may occur simultaneously. A statistical 
number of all the instantaneous bubbles (in single and multiple form) is designated as 
n OTr . The size of the bubbles appearing on this surface will vary from site to site, and 
the size of each individual bubble will vary during the growth period; nevertheless, a 
statistical bubble, which represents an average size of a typical bubble growing on the 
surface, can be envisioned. Such statistical information can be obtained by tediously 
examining high-speed motion pictures of boiling on a surface. An instantaneous statis- 
tical picture of the fraction of the surface covered by vapor can be determined. 

For the model presented herein, the projected surface area underneath the bubbles 
is of primary interest. The general bubble shape is assumed to be that of a truncated 
sphere; thus, the contact surface area between the bubble and the surface may vary 
anywhere from a circle of the bubble radius r^ (hemispherical bubble) to a circle of 
only a fraction of the bubble radius. It is assumed that the large multiple bubbles can 
also be represented by this geometry. 

A profile view of three statistical bubbles (n ay = 3) on the surface is presented in 
figure 5(b). The average contact area of the bubbles can be related to a contact angle /3, 
defined in the figure. This contact angle is often associated with the liquid contact angle 
measured on the meniscus of capillary tubes. Although the bubble contact angle and the 
liquid contact angle may be related, they are not identical. Measurements of bubble 
contact angles appear in the literature. References 34 and 35 show that the contact angle 
of a bubble varies during the growth time. For the model presented herein, a value of 
(3 averaged over the growth period of a bubble is used, and isolated bubble values for /3 
are applied to the multiple bubbles. 

Figure 5(a) shows phantom drawings of areas that will eventually be populated with 
bubbles. These areas were drawn equivalent to the projected bubble area A^. The 
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thermal layer is being conditioned during a waiting period for the generation of new 
bubbles. The ratio of the statistical area of growing bubbles to the area in preparation 
is a function of the ratio of the growth to the waiting periods. If the waiting period is 
longer than the growth period, more surface area of preparation than area of growing 
bubbles is present instantaneously. The relative distribution of the bubble and prepara- 
tion areas is reversed when the growth period is greater than the waiting period. (In 
fig. 5, the waiting period is twice the growth period. ) At higher heat fluxes, more sites 
would become active. 

Nonboiling area . - The area of the surface (fig. 5(a)) that is outside the solid line, 
or phantom areas shown, is considered a nonboiling area. This area of nonboiling may 
be appreciably large at low values of the average boiling heat flux and will probably dis- 
appear at heat fluxes corresponding to the critical or burnout value. 

Heat-Transfer Mechanisms 

In the previous section, a boiling heat -transfer surface was divided into (1) a 
projected surface underneath a growing bubble, (2) a surface where the thermal layer is 
being prepared for ebullition, and (3) a surface where no boiling activity takes place. 

The heat -transfer mechanisms in process over each of these areas are discussed, which 
necessitates subdividing these areas further (shown in fig. 5(a)) for the specific mecha- 
nisms involved. A pictorial description of the various heat -transfer mechanisms that 
make up the overall heat transfer is presented in figure 6. Essentially, this is figure 5 
modified to illustrate where the various heat -transfer mechanisms are in effect. 

Bubble-Area Mechanisms 

As shown in figure 6, the form of the bubble over this area is assumed to be a 
truncated sphere. The surface area in contact with the bubble is assumed to be wetted, 
with the exception of a possible dry spot that develops during the bubble lifetime. The 
growth of the dry spot is discussed in reference 33. The behavior of the dry spot is dis- 
cussed more fully in the appendix. Limited evidence indicates that the dry fraction of 
the bubble contact area remains constant over a broad span of the nucleate heat flux. As 
the heat flux approaches the critical, or burnout, condition, the fraction of dry area to 
total area asymptotically approaches 1 (see fig. 7). The heat -transfer mechanism above 
the dry spot is assumed to be gaseous convection, and the heat -transfer contribution is 
negligibly small. 

The wetted contact area is considered to be a microlayer (see refs. 32 and 33) that 
is the evaporative surface for the supply of vapor into the bubble. Kinetic theory from 
an evaporative surface is employed to compute the local heat transfer. 
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Questions arise regarding the legitimacy of applying the kinetic theory expression to 
the evaporative process on a microlayer surface. Generally, the picture of static condi- 
tions and thermal equilibrium at the vapor-liquid interface pertains to the kinetic equa- 
tion, but this is not the case for the microlayer in a bubble environment. Presumably, 
the microlayer is being supplied fresh liquid during the growth period, and the vapor in- 
side the bubble is highly agitated, which results in rather dynamic, nonequilibrium con- 
ditions at the liquid interface. However, the velocities of the molecular species partici- 
pating in the phase change are several orders greater than the macroscopic fluid veloci- 
ties of the bubble interior or of the microlayer. Thus, the application of equation (2c) 
to the microlayer may be a good approximation of the evaporative process occurring 
there. 

The wetted, or evaporative, area (from fig. 6) can be described as 

A evap = “av^b' 1 ' D P si ” 2 13 < 3 > 

where D r is the dry faction of the bubble contact surface. When equation (2c) is used, 
the heat transfer from the wetted surface is 


Qevap = aX( Pl ~ M 


g c M 


1/2 


,2ttRT 


sat) 


evap 


( 4 ) 


For isolated bubbles, when the bubble contact angle is less than 90°, an annulus is 
present underneath the projected bubble area that is not a part of the wetted area. This 
part of the projected bubble area and an annular area outside the projected bubble area 
(see fig. 6) are assumed to be areas of enhanced free convection. As observed in refer- 
ence 10, this enhanced region outside the bubble appeared to extend one bubble diameter 
away from the center of the bubble. The enhanced free convection is greater than the 
turbulent -free convection computed for a horizontal flat plate. Thus, a factor K is 
introduced to enhance the computed free -convection value, which applies to the convec- 
tion region under the bubble (if it exists) and the area of influence. For ordinary turbu- 
lent free convection, the equation applicable to horizontal surfaces is 


q 


= C(GrPr) 1 ^ 3 - (T 
6 w 



(5a) 


where 

C empirical coefficient for ordinary free convection 
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Gr Grashof number, V6^(T w - T^g p^/p.^ 
g local acceleration of gravity 

6 characteristic thickness 

Pr Prandtl number 

Noting that the characteristic thickness 6 cancels out and collecting similar terms 
yield 



Grouping the coefficient C and the fluid property term produces the simple empiri- 
cal equation 

Vmv = «< T w - T b> 4/3 < 5c > 


where 



For turbulent free convection on horizontal surfaces, data are available, and the 
value of 6 can be determined. Thus, the heat -transfer contribution for the surface 
involved in enhanced convection is 

Q = K( lconv n av 7Tr b (4 " sin2 ® < 6 ) 

The selection of a value for the factor K is a problem. The average effects of 
superposing bubbling on a free -convection process were considered in reference 18. 
Electrolytically generated bubbles produced enhancement ratios that ranged from slightly 
greater than 1 to as high as 4. The latter value occurred at the lower temperature 
differences, and, consequently, may reflect changes in the free -convection mechanism 
from a laminar to a turbulent free convection. With boiling present, a laminar-type free 
convection does not seem likely; thus, the enhancement ascribed to bubbling would be 
more comparable to what the authors of reference 18 saw at higher temperature differ- 
ences, where the maximum enhancement ratio was approximately 2. As an average 
representative value of enhancement, 1. 5 would be a reasonable approximation. Con- 
sequently, a constant factor of 1. 5 is the multiplier K of equation (6) to account for 
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the enhanced convection ascribed to bubble-growth effects. This value is assumed to 
apply to all fluids, although it should be a variable and account for the different rates 
of bubble growth, the size of bubbles, and the density of bubble sites - to mention a few 
influences. Because no direct way of accounting for these variations is apparent, the 
enhancement factor is considered to be a constant in the model proposed herein. 


Area of Preparation Mechanism 

Shown in figure 6 are areas where the thermal layer was making a transition to a 

condition that would produce ebullition. As discussed in reference 22, this period of 

transition, commonly known as the waiting period, can be approximated as a transient 

conduction process into a semiinfinite fluid slab. The equation for transient conduction 

is well known. The length of time involved in the transient conduction is taken as the 

waiting period f , and the contribution from the area of preparation is 
w 

Qcond = (Tw ' Tb>Aprep <7> 

2 

The area of preparation A__-_ is n^fr/r) nrf. It was mentioned earlier that, 

prep dv w gr o 

if all the sites could be typified by one identical type, the ratio of the preparation areas 
to the growth areas would be proportional to the ratio of the waiting to the growth periods. 
At the higher nucleate heat fluxes, a situation will arise in which the areas of prepara- 
tion will overlap the area of evaporation. This situation is discussed in detail in the 
appendix. The temperature difference for the waiting period was assumed to be between 
that of the wall temperature and the bulk conditions in the liquid. 


Nonboiling Surface Mechanism 

The nonboiling surface constitutes the area remaining on the surface that is not 
involved directly in any aspect of the ebullition process. The mechanism of heat transfer 
is turbulent free convection. Equation 5(c) describes this process. The area involved 
in the free -convection process is 
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A - a - A -A -A -A 

conv h evap enh conv prep d 


jj. 

= A h " n av ffr b^ ' Dr ) sin2 0 “ n av 7rr b( 4 ~ sin2 ^) " n ^ v nr h — " n av 7rr b Dr sin2 ^ ( 8 ) 


gr 


Q 


conv 


= ( ’< T W - T b> 


4/3 


A h - “av^bl 4 + T W 

gry 


(9) 


General Equation for Model 


An equation for all the heat -transfer mechanisms in process at any instant over the 
entire surface is simply the sum of all the contributions for each of the subareas. When 
this equation is written as a heat flux a , it is expressed as follows: 


Q _ 


= < W = 


(1 - Dr)sin 2 /3 


L A h 


Mp* - Pr> 


g^ 

2ttRT 


, 1 / 2 * 


+ 


\v n l 


(4 - sin'p) 


K 0 ( T w - T b ) 


4/3 


r n nr? r 1 
av b w 

A h T gr/ 


2 k- 




< T w - T b> 


II CUT, 


W 


1 _ av b [4 + .'W 


h 


gr/J 


e(T w - T b ) 4/3 (10a) 


In the actual calculation procedure, it was convenient to use equation (10a) in the 
following form: 
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lav = t 4 * 1 * Dr > sin2fJ ]^vap + [ 4 < 4 - sl " 2 «] K( l CO nv 


+ # 


w 


^cond 


1 - # 4 + 


gr 



conv 


(10b) 


The application of equation (10a) to compute an average overall heat flux requires 
considerable experimental input information about the particular boiling condition in 
question. Statistical information on the amount of surface area covered by bubbles and 
the relative magnitudes of the growth and waiting periods is required. Surface and bulk 
temperature data are also needed. Thus, the model proposed herein cannot be considered 
to be a prediction technique, which would be capable of predicting the fraction of surface 
area covered by bubbles and the relative magnitudes of the growth and waiting periods. 
Predicting this area is a difficult task that requires more detailed information than is 
presently known about the formation and growth of bubbles. Further studies are needed 
on the evaporation, condensing, and convection processes associated with establishing 
bubble size. Vapor appearing in a bubble is speculated to be the instantaneous surplus of 
vapor not condensed over the bubble cap. If this is true, bubble size could be predicted 
from information about the efficiency of evaporation and condensing rates over the 
bubble surface. 

Equation (10a) was developed as a means to compute the relative contributions of 
several heat -transfer mechanisms considered to be important in the overall nucleate - 
boiling heat transfer. Such a calculation will lead to an assessment of the relative 
importance of the heat-transfer mechanisms for a known nucleate -boiling condition. 

The appendix contains calculations made by employing equation (10a) for nucleate - 
boiling data in which temperature data and statistical information on the fraction of the 
surface covered by single or multiple bubbles were available. Water and methanol 
were the fluids used in the comparison between the computed and measured heat fluxes. 

References 11 and 36 were used to obtain heat flux and temperature data for water. 
Reference 37 was used to obtain heat flux -temperature information for methanol. 
Reference 38 provided information on the average bubble population (single and multiple 
bubbles) and on the average surface area covered by these bubbles (for both water and 
methanol). These data were obtained from a careful examination of high-speed motion 
pictures and represent a statistical analysis of many data samples. Reference 38 
presents information on the average heat-flux and bulk-fluid conditions but gives no 
surface -temperature data. Reference 39 contained near-burnout data for water. 
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Factors Not Considered in Model 


Since the proposed model involves transient processes near the liquid-solid inter- 
face, the smoothness of the solid surface, its thickness, and the thermal diffusivity of 
the material will be influential. Such effects are not included in the equations for the 
model presented herein. 

Surface smoothness . - Numerous observations have been made which show that 
nucleation is strongly related to the surface smoothness. Hsu (ref. 22) demonstrated 
that the thermal conditions (thermal -layer profile and thickness) for a nucleation 
criterion can be related to the dimensions of a scratch or pit. However, no one has 
been able to present a general correlation for bubble population, frequency, and surface 
distribution as a function of surface roughness. Thus, an accounting for the surface 
cannot be included in the model presented herein. 

Surface thickness and geometry . - Experimental boiling data that cover a variety of 
heater geometries, including wires, thin ribbons, rods, tubes and thick blocks are 
presented in the literature. These heater geometries can influence the surface -tempera- 
ture distributions and the transport of heat into the fluid. For example, a thin metal 
ribbon would differ significantly from a thick block of the same material in that the block 
would be apt to provide a constant surface temperature and a variable instantaneous heat 
flux, whereas the thin ribbon would approximate a constant heat flux with a variable 
surface temperature. 


Thermal Diffusivity of Metal Surface 

Thermal balance across the solid -liquid interface is influenced by both the transient 
and steady modes of transport within the metal heater. The transient effect is discussed 
in reference 33, where it is shown that the transient thermal transmissivity of a finite 
plate is a function of k/ yfa. This parameter was computed to range from 0. 03 calories 
per cubic centimeter per °C for Inconel and glass to 0. 09 calorie per cubic centimeter 
per °C for copper. Thus, copper will maintain a more constant surface temperature 
than glass or Inconel. The heat that is transported into the fluid at the bubble area is 
influenced by the thermal diffusion of the heater (ref. 6), Consequently, the diffusion 
properties of the heater surface are important during transient conduction into the 
bubble thermal layer. 
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TABLE m. - BASIC EXPERIMENTAL INFORMATION FROM REFERENCE 36 
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RESULTS AND DISCUSSION 


The model of the nucleate -boiling mechanism described herein was utilized in a cal- 
culation of the various components of heat transfer that were assumed to be involved in 
the overall nucleate -boiling process. The input data employed in the calculations are 
given in table in. As mentioned in the section Description of Model, the assumptions 
and approximations related to the use of these data are part of the appendix. 

The results of the calculation are presented in table IV. The contribution of each 
heat -transfer mechanism is tabulated along with the overall average heat flux. The 
final column in the table presents the experimental average heat flux from reference 38. 
Note that the tabulations include computed time-averaged values and values where both 
time and heat -transfer surface weighting are involved. The dashes in the table signify 
that area overlapping of mechanisms associated with the ebullition process appropriated 
the entire surface and left no areas for convection. The procedure used for this situation 
in the calculation is discussed in the appendix. 

At first glance, the comparison of experiment with analysis in table in may seem 
limited in terms of the number of comparisons. However, it should be realized that the 
source of information on the bubble areas, the average growth and waiting periods, the 
number of bubbles, etc. used in the comparison represents an exhaustive statistical 
study of hundreds of bubbles (ref. 38). Realistically, these six runs should be con- 
sidered as representative of scores of nucleate -boiling data points. 

The last two columns of table IV do show remarkable agreement in the computed 
overall heat flux and the experimental values. It should be noted that the calculations 
for water represent relatively low heat fluxes with respect to critical or burnout heat 
flux. Information on the surface area involved in the merging bubbles for water could 
not be discerned easily; consequently, reference 38 did not contain information relevent 
to near-critical heat fluxes of water. Reference 39, however, did present instantaneous 
measurements of the bubble contact area over a range of fluxes including the critical 
value. No data were presented on the number of sites, the merging of bubbles, or the 
bubble frequency. The ratio of bubble growth to waiting time was computed by assuming 
that 


T w s A prep ^ 

T gr A evap 

subject to the constraint that (A__ a _ + A„ ro „)/A. < 1 (see appendix for further discussion). 

prep Gvdp n 

These near -burnout comparisons for water (based on ref. 39) are tabulated separately 
in table V for two heat -flux levels at saturated and subcooled conditions. Again, the 
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computed results agree reasonably well with the experimental measurement. For some 
unknown reason, the agreement is better for the subcooled cases. 

The stated objective of this report was to present a balanced observation of the 
relative contributions of several mechanisms that are thought to make up the overall 
nucleate -boiling heat flux over a range of heat fluxes. In figure 8, the relative contribu- 
tion of each heat-flux component is plotted as a function of the percent of critical heat 
flux. An examination of the figures leads to the following observations: 

(1) For both water and methanol, no single mechanism dominates over the entire 
boiling heat -flux range. It is true, however, that evaporation becomes the chief con- 
tributor for heat fluxes that are greater than 20 percent of the critical value. The well- 
known hydrodynamic -instability, critical-heat-flux models of Zuber and Kutateladze 
(see ref. 40) indicate that evaporation is the dominant heat-transfer mechanism near the 
critical heat flux. In their model, the critical heat flux develops when the hydrodynamic 
instability limits the amount of heat -transfer surface wetted by the liquid phase. In the 
analytical expression, the critical flux is directly proportional to the heat of vaporization, 

(2) The relative contribution of the evaporation mechanism is higher with water 
than with methanol. Since water does have a much higher heat of vaporization than 
methyl alcohol, this trend would be expected. 

(3) Transient conduction is an important contributor to the total heat flux over a 
broad range of conditions. The cyclic removal of heat from the surface into the liquid 
sublayer by conductions appears to be a very efficient process compared with steady- 
state free convection. The repeated destruction of the thermal layer that ensues the 
departure of a bubble enables this transient -conduction process to be repeated at a 
bubble site. 

The overall model presented herein emphasizes that nucleate boiling should be 
thought of primarily as a series of transient heat -transfer processes over distinct 
regions of the surface. The computation of an average heat flux requires knowledge of 
the time and spacial averaging of the heat -transfer processes. 

Considerable experimental information is still needed to give more comprehensive 
information on aspects of the nucleate -boiling mechanisms over a range of heat flux and 
subcooling. The limited amount of information on the bubble contact area needs to be 
expanded. The behavior of the microlayer and the dry spot should be studied at higher 
heat fluxes. The distribution of condensing and evaporating heat-transfer coefficients 
over a bubble surface requires study. The diffusion of the thermal energy into the fluid 
bulk is another subject that deserves some attention for better overall understanding. 

The model discussed herein treats the heat transport across a fluid resistance at the 
interface between the surface and the fluid. The means for absorbing that thermal 
energy in the bulk is not considered. In two-phase flow processes, the mechanism of 
thermal-energy absorption in the bulk is particularly important. The^void distribution 
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across the channel and the enthalpy gradients in both the radial and axial directions are 
related to it. 

Mention was made earlier of the effect of the surface condition and the surface 
material on the overall boiling process. The limited experimental information on this 
subject suggests that this topic be further pursued. 

In this discussion, several obvious deficiencies in experimental information on 
nucleate boiling were cited. The list is by no means complete. Even a partial listing 
gives emphasis to needed research in nucleate boiling. 


SUMMARY OF RESULTS 

1. An overall model of nucleate boiling, based on recent experimental observations 
by several investigators, is presented. The model incorporates convection, transient 
conduction, and evaporation - heat -transfer mechanisms. The contribution of each 
mechanism is weighted according to spacial and time averaging. 

2. Statistical experimental data for water and methanol were used as input to the 
heat flux equations of the model. The resulting computed overall heat fluxes agreed 
well with experimental results over a wide range of heat flux, 

3. No one single mechanism dominates over the entire range of overall nucleate - 
boiling heat flux. Evaporation is the major contributor when the heat flux is greater 
than 20 percent of the critical or burnout value. Transient conduction is next in 
importance . 

4. The use of the analytical model emphasizes the need for more basic information 
on aspects of nucleate boiling over a range of heat flux and subcooling. Studies on bubble 
contact area, evaporative -condensing processes in bubbles, and thermal diffusion into 
the liquid buik would make possible predictions of the input information needed for the 
model. Such generalized information may transform the model from just an analytical 
description of nucleate boiling into a prediction technique. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, January 5, 1967, 

> 5 129-01-09-04-22. > | 
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APPENDIX - CALCULATION OF BOILING HEAT FLUX 
FROM POSTULATED MODEL 


It was pointed out previously that experimental information is needed to compute the 
average boiling heat flux from equation (10a), The model presented herein will not 
predict the surface conditions. Experimental information from references 11 and 36 
to 38 is employed to compute heat flux for water and methanol. This input information 
is given in table in. In addition to the assumptions and approximations included in the 
model description, certain other assumptions and approximations are employed in the 
calculation that are explained in this appendix. Some repetition of the assumptions 
mentioned previously are necessary. A sample calculation appears at the end of this 
appendix. The various heat -transfer mechanisms associated with the division of areas 
shown in figure 6 are discussed separately. 


Evaporative Area Under Bubble 

Dry area . - As depicted in figure 6, the evaporative area is the base of the truncated 
spherical bubble minus the dry -spot area. The fraction of the bubble contact surface that 
was dry was assumed to be a function of the magnitude of the average boiling heat flux. 
Figure 7 presents a general curve in which the dry fraction of the contact area was one 
coordinate and the fraction of burnout heat flux was the other. Information from refer- 
ences 33, 37, and 39 permits the construction of a general curve of this type. Fortu- 
nately, the type of heating surface employed in each of these investigations was similar 
(conductive glass). Surface tension has not been considered in developing figure 7. 

The circles in figure 7 represent data obtained from reference 33 for heat flux, 
which are approximately 2, 3, and 4 percent of the critical value (the critical heat flux 
for water is taken as 5xl0 5 Btu/(hr)(ft 2 ) (1, 57xl0 6 W/m 2 ). Figure 16(b) of reference 39 
shows an instantaneous view of dry spots at a heat flux which is at 20 percent of burnout. 
The dry area is similar in shape and magnitude to that of the lower heat fluxes in refer- 
ence 33. The square symbol in figure 7 is an estimate from figure 16(b) of reference 
39. The authors of reference 37 observed that the dry-spot area fraction remained 
substantially constant up to 70 percent of the critical flux. Consequently, the portion of 
the curve in figure 7 for heat flux, which varies from 4 to 70 percent of the critical flux, 
is drawn with the percent dry spot as a constant. 

The onset of burnout was described in references 37 and 39 as a highly transient 
state at the local burnout point where film boiling rapidly replaces multiple -bubble 
nucleate boiling. Thus, above 70 percent of the burnout flux, the fraction of dry area 
for a constant heat-flux surface grows asymptotically toward 1. 
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Contact angle . - Reference 31 presents contact -angle data for steam bubbles in 
water. The average contact angle was 30°. No published data were found for methanol; 
however, it is expected that the contact angle would be smaller because the surface 
tension of methanol is less than that of water. An estimated value of 20° was used in 
the calculations. High-speed motion pictures of ethanol (private communication with 
R. Siegel and E. G. Keshock) indicated that the contact area of the bubble could be 
approximated by a 20° contact angle. 

The experimental information on the contact angle was obtained from isolated 
bubbles. It was assumed that large multiple bubbles could be considered to be 
agglomerates of single bubbles and that the same contact -angle information which was 
applied to the multiple bubbles could be applied to the single bubbles. 

Pressure difference in evaporation . - Equation (4) contains a pressure-difference 
term p^ - p r , where p^ was assumed to represent a saturation pressure correspond- 
ing to an average temperature of the microlayer during the bubble growth 



where T is the saturation temperature at bulk conditions of liquid. This average 
temperature was chosen because figure 3 shows that the surface temperature under a 
bubble undergoes an almost linear drop during the early stage of the growth period when 
evaporation is significant. An arithmetic average is a good approximation to an average 
surface temperature, and the microlayer temperature must be related to the wall 
temperature. 

The inherent sensitivity of the saturation -pres sure value to the saturation temper- 
ature influences greatly the term p^ - p r in equation (4). The selection of p^ is an 
empiricism that can be judged on the reasonableness of the estimate of Q_„ 0 « obtained 
near the burnout condition. The numerical results (table V) indicate that the previously 
described method for selecting p^ is reasonable. The evaporative coefficient a in 
equation (4) was obtained from the literature (see table I, p. 9). 

Surface of Boiling Preparation 

This portion of the surface is assumed to be undergoing transient conduction while 

an area around a site is prepared for the birth of a bubble. The area in preparation is 

o 

assumed to be equal to the statistical bubble area nr^. 

Equation (7) requires experimental information on the waiting period. The average 
waiting period can be computed form the statistical data of a boiling experiment in the 
following fashion. The ratio of the number of instantaneous bubbles n ay to the number 
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of bubble sites N of a surface (see table HI, p. 21): 


gr 


n, 


av 


t + T 
w gr 


N 


or 


w 


T + T 

w gr 


= 1 


n 


av 

N 


( 12 ) 


Obtaining a numerical value for t w requires that a value be assigned to Tg r . It is well 
known that the growth time for individual bubbles is a function of heat flux and subcooling. 
In figure 12 of reference 37, the dependence of growth time on heat flux is shown for 
carbon tetrachloride. The data encompass the entire nucleate -boiling regime from 
incipience to burnout. Herein, this dependency of growth time on heat flux for carbon 
tetrachloride (ref. 37) was generalized, and a curve of a normalized growth time, 

TgrAgr re j was plotted as a function of the percent of the critical heat flux in figure 9. 
For the fluids considered, ref was assumed to be 10”^ second. This value is 
representative of much experimental data for water, but it is probably too large for 
methanol. The shaded region in figure 9 reflects the uncertainty in the values obtained 
for heat flux and bubble -growth times in the vicinity of the incipience point. No attempt 
was made to correct the bubble -growth times for subcooling effects. 

For ease of calculation, the thermal properties inserted in equation (7) were taken 
at bulk conditions. Perhaps a more accurate selection would be to use film properties, 
but the result would not differ appreciably. 

As pointed out previously in this report, the area of preparation A is assumed 

prep 

to be related to the ratio T w / T g r and to the area of evaporation A gva p. At the higher 
heat -flux portion of the boiling curve, the areas of preparation and evaporation of 
multiple bubbles probably overlap. It was assumed in this model that no such overlap 
occurs. Equation (11) (and its restraint) is valid over the entire range of the nucleate- 
boiling heat flux. An inconsistency between the ratio of evaporative to preparation 
areas and the ratio of growth to waiting times appears in run 5 (table III) and is attributed 
to overlapping. For this near burnout run, r w / T g r is much greater than ^- eva p/Ap re p 
or cp/(\ - <&) and thus violates the restraint of equation (11) [$ + ( T w / T g r )(®) < 1)] • 

For the burnout and near -burnout runs of reference 39, the restraint of equation (11) was 
not violated. No statistical information on bubble sites and population was available; 
therefore, T w / T g r was estimated from information on the instantaneous areas of 
evaporation and from equation (11) with the restraint applied. 
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Nonboiling Surface 


The estimation of the turbulent free -convection heat flux was discussed in the 
previous section of this appendix. The area undergoing turbulent free convection is 
assumed to be that surface which is not participating in some aspect of the ebullition 
process or is not directly linked to the bubble growth, as is the case for the area of 
influence (enhanced turbulent convection). 

For some of the computations with methanol, no surface area was left for the free- 
convection component because the entire surface was engaged in the boiling mechanisms. 
In fact, for run 5, so much of the surface was engaged in supporting ebullition that no 
enhanced or free -convection areas were available. What was done in these cases was 
to start the computation with the evaporative area. The residue of area left to other 
mechanisms was taken in the order of transient conduction first and then enhanced con- 
vection. 

Justification for the assumption of the overcrowding of a surface by the boiling 
mechanisms can be seen in reference 35, wherein it was noted that the spacing between 
active sites was reduced rapidly as the heat flux was increased. 


Sample Calculation 

The following sample calculation presents the estimates of each contribution to the 
overall heat flux for nucleate boiling in water. Similar calculations were made for 
methanol. 

The following input data were inserted in the heat flux equations of the model for 
water: 

(1) Properties of water (including evaporation coefficient): 

T sat - 212 ° F 
P = 1 atm 

a = 6. 5xl0 -3 ft 2 /hr 
k = 0. 39 Btu/ (hr) (ft) (°R) 

X = 970 Btu/lb 


(2) Statistical boiling information: 


a = 0.04 
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From references 11 and 38, T - T,= 27° F and T = 225° F. Since 

W SJ W 

T g r / ( T g r + T W ' = 16 ’ and T gr = T w + T gr = °* 06 second 5 thus, = 0. 05 second 

Evaporation contribution (first term of eq. (10)): 

/ ScM Y /2 

Wp = ^-Pr)(j— J 

where is the saturation pressure at 

T + T j. 

T = — = 218. 5° F 

2 

p^ = 16. 8 psia 

q pva = (0. 04)(970)(16. 8 - 14. 7)144^— 18 X 32 ' 2 3600 

P \2ttX 1544 x 672/ 


= 400 000 Btu/(hr)(ft 2 ) 


Conduction contribution (third term of eq. (10)): 


^cond = 2k < T w 



= 2x0. 39(27) 


' 14x10' 



= 34 500 Btu/(hr)(ft 2 ) 


Convection contribution: 

Rather than computing the free -convection contribution, recourse is made to refer- 
ence 36 where indications of the free -convection heat flux were measured. 

%nh conv = 11 000 x 1 - 5 = 16 500 Btu/(hr)(ft 2 ) (second term of eq. (10)): 

The multiplier 1. 5 in q 0nh cony is the enhancement factor. 
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q free conv = 11 000 Btu/(hr)(ft 2 ) (fourth term of eq. (10)): 
Area-weighted contribution: 


O 

^evap 



q evap 


A 

evap 

A h 


= 400 000 [$(1 - Dr)sin 2 ^] 


For 


/3 = 30° 


and 


Dr = 0.15 


Q 

evap = (400 000)(0. 0173) 
A h 


= 7000 


^cond 


q cond ' 


cond 



= 45 500 (0.436) 


= 14 900 


(a) 


(b) 
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Qenh conv _ A enh conv 

» ^enh conv A 

A h A h 

= 16 500 [#(4 - sin 2 0)] 

= 16 500 (0. 312) = 5100 

Q A 

^conv _ conv 

A q conv A 

n A h 

■ 11 000 & - (A evap + A cond + A enh co n v>] 
= 11 000 (0.235) 

= 2600 


The summation of all the contributions is 


Z -Q- = 29 600 Btu/ (hr) (ft 2 ) 
A h 


(c) 


(d) 


(d) 
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(a) Bubble agitation. 


- Hot liquid 
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(b) Vapor-liquid exchange. 
Condensing 



^Microlayer 
<c) Evaporation. 

Figure 1. - Schematic diagram of nucleate boiling models. 
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Thermal 


Thermal layer 


Heater surface 


Established thermal layer 


Bubble 


Bubble growth period 




Destruction of thermal layer 
figure .2, - Principal stages in life of thermal layer. 









Temperature, 



Figure 3. - Temperature and heat flux beneath growing bubble. Saturation temperature, 95° F ; liquid bulk temperature, 
84° F; liquid, water; reduced pressure. 
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\ by bubbles — 


Nucleation sites - 



(a) Plan view of surface. 


Areas of preparation 



Figure 5. - Instantaneous representation of nucleate-boiling surface. (The 
waiting period is equal to twice the growth period.) 
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(b) Profile view. 

Figure 6. - Distribution of heat-transfer mechanisms over nucleate-boiling surface. 
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mass), D r 


Figure 7. - Fraction of dry surface under bubble as function 
of dimensionless heat flux. 
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Free convection 
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(a) Water, subcooled. 
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Percent of critical flux 


(b) Methanol, subcooled. 

Figure 8. - Relative contributions of heat-transfer 
mechanisms over range of average heat flux. 
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Percent of critical flux 

Figure 9. - Growth time of bubbles as function of 
percent of critical heat flux. (Based on fig. 12 
from ref. 37.) 
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Motion picture C-252 is available on loan. Requests will be filled in the order 
received. You will be notified of the approximate date scheduled. 

The film (16 mm, 7 min, color, sound) shows actual high-speed motion pictures of 
evaporative microlayers underneath growing bubbles. Reflective images and Newton 
rings of the microlayer are presented as evidence of the existence of the microlayer. 

Motion picture C-252 is available on request to 

Chief, Technical Information Division (5-5) 

National Aeronautics and Space Administration 
Lewis Research Center 
21000 Brookpark Road 
Cleveland, Ohio 44135 



CUT 

Please send, on loan, 

Date 

copy of Motion Picture C-252 to 

Name of organization 


Street number 



City and State Zip code 

Attention: Mr. 
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tr The aeronautical and space activities of the United States shall he 
conducted so as to contribute . . . to the expansion of human knowl- 
edge of phenomena in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof” 

—National Aeronautics and Space Act of 1958 
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